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ABSTRACT 

We present dust continuum observations of the protoplanetary disk surrounding the pre-main sequence 
star AS 209, spanning more than an order of magnitude in wavelength from 0.88 to 9.8 mm. The 
disk was observed with sub-arcsecond angular resolution (0.2" — 0.5") to investigate radial variations 
in its dust properties. At longer wavelengths, the disk emission structure is notably more compact, 
providing model-independent evidence for changes in the grain properties across the disk. We find that 
physical models which reproduce the disk emission require a radial dependence of the dust opacity n v . 
Assuming that the observed wavelength-dependent structure can be attributed to radial variations 
in the dust opacity spectral index (/?), we find that f3(R) increases from f3 < 0.5 at ~ 20 AU to 
P > 1.5 for R > 80 AU, inconsistent with a constant value of f3 across the disk (at the 10a level). 
Furthermore, if radial variations of n v are caused by particle growth, we find that the maximum 
size of the particle-size distribution (a max ) increases from sub- millimeter-sized grains in the outer 
disk (R > 70 AU) to millimeter and centimeter-sized grains in the inner disk regions (R < 70 AU). 
We compare our observational constraint on a max (R) with predictions from physical models of dust 
evolution in proto-planetary disks. For the dust composition and particle-size distribution investigated 
here, our observational constraints on a max (R) are consistent with models where the maximum grain 
size is limited by radial drift. 

Subject headings: protoplanetary disks — stars: individual (AS 209) 



1. INTRODUCTION 

The growth of sub-micron-sized dust grains into 
millimeter and centimeter-sized particles is a funda- 
mental component of the planet formation process 
(Beckwith et al. 2000; Natta et al. 2007). Grain growth 
directly affects the optical properties of dust parti- 
cles (Henning & Stognienko 1996; D'Alessio et al. 2001), 
with composition and temperature having a lesser im- 
pact (Henning & Mutschke 1997; Semenov et al. 2003; 
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Boudet et al. 2005). At long wavelengths, the slope 
a of the spectral energy distribution (SED, vS v oc 
v aJrl ) can be related to the spectral index of the dust 
opacity f3 [k v oc in particular f3 = a — 2 for 

optically thin warm dust (Beckwith & Sargent 1991; 
Miyake & Nakagawa 1993). The millimeter dust opac- 
ity slope for small interstellar medium (ISM) dust grains 
is PlSM ~ 1-7 (Li & Draine 2001). A smaller value of f3 
is measured if larger millimeter-sized grains are present 
(D'Alessio et al. 2001; Draine 2006). 

Observations of dust in disks from sub-mm to 
cm wavelengths (Testi et al. 2001; Calvet et al. 2002; 
Testi et al. 2003; Natta & Testi 2004; Wilner et al. 2005; 
Rodmann et al. 2006; Lommen et al. 2009; Ricci et al. 
2010a,b) have provided strong evidence for grain growth 
in disks. These results, obtained mainly from spa- 
tially unresolved observations, assumed a constant dust 
opacity (i.e. constant j3) throughout the disk. How- 
ever, theoretical models of dust transport, fragmenta- 
tion, and size evolution in proto-planetary disks pre- 
dict that the average size of grains varies with the dis- 
tance to the central star (hereafter referred to as ra- 
dius, Dullemond & Dominik 2005; Birnstiel et al. 2010). 
Hence, radial variations in the dust opacity slope are ex- 
pected. 

The increased sensitivity and angular resolution of 
current interferometers has enabled the first studies of 
radial variations of grain growth within circumstellar 
disks. Observations of disks in Taurus at 1.3 and 
3 mm with CARMA 15 (Isella et al. 2010) and PdBI 16 
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(Guilloteau et al. 2011), have constrained f3 as a function 
of radius. However, these results are only sensitive to ra- 
dial variations larger than A/3 ~ 0.6 — 0.7 at 3<r in the 
best cases, limited by the small wavelength separation 
between the two bands being observed. Increased wave- 
length coverage has been obtained by adding cm-wave 
observations (Banzatti et al. 2011), but the low signal- 
to-noise ratio, before the Karl G. Jansky Very Large 
Array (VLA) upgrade, made such studies challenging. 
To significantly improve these constraints, and effectively 
compare them to theoretical predictions of grain growth, 
increased wavelength coverage and enhanced sensitivity 
are required. 

This letter presents the first results of a program 
that combines spatially resolved observations from sub- 
millimeter to centimeter wavelengths, to constrain ra- 
dial variations of grain growth within circumstellar disks. 
The pre-main sequence star AS 209, also known as VI 121 
Oph, is a classical K5 T-Tauri star with high accre- 
tion rate (1.3 x 10" 7 M yr _1 , Johns-Krull et al. 2000). 
AS 209 is isolated from the main p-Ophiuchus cloud 
(a = 16 /l 49 m 15.3 s ,(5 = -14 d 22 m 08.7 s ) at a distance of 
131 ±50pc (van Leeuwen 2007). Resolved 0.88 mm con- 
tinuum observations reveal a smooth radial density struc- 
ture down to 0.3" ~ 40 AU (Andrews et al. 2009, 2010). 
We obtained additional spatially resolved observations at 
A > 0.88 mm, to study grain growth inside the AS 209 
disk. 

2. AS 209 OBSERVATIONS 

Submillimeter Array (SMA) observations at 0.88 mm 
were obtained between May 2006 and June 2007, in C, 
E, and V configurations (~ 0.3" resolution), with total 
continuum bandwidth of 4 GHz. Its integrated flux den- 
sity is 577 ± 60 mJy at 0.88 mm. Andrews et al. (2009) 
presented a detailed discussion of the observations and 
data reduction. 

CARMA observations at 2.8 mm were obtained be- 
tween December 2009 and March 2010, in C, B, and A 
configurations spanning baselines of 30-1800 m (~ 0.5" 
resolution), with total continuum bandwidth of 8 GHz. 
Its integrated flux density is 40 ± 6 mJy at 2.8 mm. At- 
mospheric conditions were good, with zenith opacities 
r 230 GHz < 0.15, except for B configuration observations 
where T230 GHz °° 0-4- Data were calibrated using the 
Multichannel Image Reconstruction, Image Analysis and 
Display (MIRIAD) software package (Sault et al. 1995). 

VLA observations at 8.0 and 9.8 mm (Ka-band) were 
obtained between Jan-May 2011 (project AC982), in 
CnB and BnA configurations, spanning baselines of 
350 m to 36.4 km (~ 0.2" resolution). Two 1 GHz base- 
bands centered at 30.5 and 37.5 GHz provided 2 GHz of 
continuum bandwidth. At 8.0 and 9.8 mm the integrated 
flux density of AS 209 is 1.2 ± 0.1 and 0.7 ± 0.1 mJy 
respectively. Data were calibrated using the Common 
Astronomy Software Applications (CAS A). 

For all datasets, quasar observations interspersed be- 
tween science target observations were used to calibrate 
the complex antenna gains. Absolute flux calibration 
was obtained by observing Neptune for CARMA, Uranus 
and Titan for SMA, and 3C286 for VLA. For each tele- 
scope, calibrated visibilities coming from different config- 
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Fig. 1. — Real and imaginary part of the bin-averaged visibili- 
ties from AS 209 as a function of spatial frequency (uv -distance). 
Disk visibilities have been normalized by the measured flux den- 
sity between wv-distances of — 40 [kX], and deprojected assum- 
ing PA= 86°, i = 38° (Andrews et al. 2009). Filled circles and 
errorbars: correlated emission observed at different wavelengths, 
continuous lines: best-fit at each wavelength. 



urations were corrected for proper motion (van Leeuwen 
2007) and compared over the spatial frequencies where 
these configurations overlap; these were found to agree 
within the absolute flux scale uncertainty (~ 10 — 15%). 
Fourier inversion of the visibilities and cleaning of the 
dirty image to produce dust emission map, were per- 
formed using CASA. 

3. OBSERVATIONAL RESULTS 

Figure 1 presents the normalized Real/Imaginary part 
of the correlated emission from AS 209 as a function 
of spatial frequency (i^-distance), for each wavelength. 
The real part of the visibility will be constant with in- 
creasing spatial frequency for an unresolved point-source. 
In contrast, for AS 209 we measure a declining profile 
with increasing uv-distance, indicative of spatially re- 
solved emission at all wavelengths. Furthermore, short- 
wavelength emission (0.88, 2.8 mm) diminishes faster 
than long- wavelength emission (8.0, 9.8 mm); therefore, 
the structure at long- wavelengths is more compact than 
at short- wavelengths. 

Before interpreting these observational results in terms 
of grain properties variations, we explore other expla- 
nations for the differing amplitude profiles: (1) see- 
ing affecting high frequency observations would have to 
be > 1.2" at 0.88 mm (xlO larger than observed see- 
ing; Andrews et al. 2009), (2) optically thick emission at 
short wavelengths: modeling of these observations (§4) 
indicates r < 1 over the scales studied. Consequently, we 
conclude that the long-wavelength emission comes from 
a smaller region than the short-wavelength emission. 

Before proceeding, we account for continuum emis- 
sion not arising from dust, such as free- free emission 
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Fig. 2. — Left: T\(R) inferred from separate modeling of our mult i- wavelength observations, assuming a constant dust opacity with 
radius. Right: Optical depth t\(R) = k\ x T,\(R) inferred from separate modeling of mult i- wavelength observations, assuming a radially 
constant k\. Colored regions: 3a confidence interval constrained by our observations, continuous line: best-fit model, dashed line on left 
panel: average temperature profile. The different ^(R) and T(R) profiles for each wavelength portray a varying dust opacity with radius, 
and because of this, none of them is the true surface density and temperature profile of the disk. 



from photo-evaporative stellar winds. Such emission 
originates near the star (within 10 AU for TW Hya, 
Pascucci et al. 2011), appearing as a point-source at the 
resolution of our observations. Indeed, in addition to the 
resolved emission detected with VLA, we detect a point- 
like component (at ^distances > 900 kA) with a flux 
density of 0.08 ± 0.02 mJy (8 mm) and 0.10 ± 0.02 mJy 
(9.8 mm), corresponding to ~ 7% and ~ 14% of the total 
emission at these wavelengths. However, VLA observa- 
tions at 5.2 cm from the Disks@EVLA collaboration 17 
(Chandler et al. in prep) place a 3a upper limit of 72/iJy, 
suggesting that the point-like emission at 8 and 9.8 mm 
may be dust emission. Nonetheless, we conservatively 
include this point-like emission as a free-free component 
in models described next. 

4. MODELING RESULTS 

Observations were analyzed using the disk emis- 
sion model described in Isella et al. (2009). The 
disk structure follows the two-layer disk approximation 
(Chiang & Goldreich 1997), defined by the boundary be- 
tween the disk surface layer (optically thin to stellar ra- 
diation) and the disk interior (optically thick to stel- 
lar emission). The disk is passively heated by AS 209 
(L* = 1.5 L ,T* = 4250 K, M* 0.9 M , Andrews et al. 
2009). Dust in the disk surface absorbs the stellar radi- 
ation and re-radiates this energy at longer wavelengths. 
Because this hot dust layer will be optically thin to its 
own thermal radiation, half of the surface layer emission 
will be radiated outwards while the rest will be radiated 
towards the disk, heating up its interior. 

For the density structure, we employ the self- 
similar solution for a viscous accretion disk 
(Lynden-Bell & Pringle 1974) 



E(R) 



x exp 



2(2 - 7 ) 

17 https://safe.nrao.edu/evla/disks/ 
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where Et corresponds to the surface density at radius R t . 
This S(i?) prescription behaves as a power- law for small 
radii, decreasing exponentially at large radii. 

We compute the dust opacity n v using Mie the- 
ory, assuming a population of compact spherical grains 
in a power-law size distribution, n(a) oc a -9 , for 
dmin < a < a max with a min = 0.01 fim. Optical 
constants for astrophysical grain components were ob- 
tained from Semenov et al. (2003), Zubko et al. (1996), 
and Warren (1984). Fractional abundances were adopted 
from Pollack et al. (1994, in volume: 8% silicates, 30% 
organics, 62% water ice). Recent revisions of the solar 
oxygen abundance (Asplund et al. 2009) would reduce 
the amount of water-ice. However, the effect on the disk 
modeling is minimal, with a difference in the disk struc- 
ture constraints well within their uncertainties as shown 
by Isella et al. (2010). 

The best-fit model to a single wavelength observation 
is found through x 2 minimization with 3 free parame- 
ters that describe the surface density (E £ ,i^t,7), since 
the disk geometry is fixed based on molecular line obser- 
vations (i = 38°, PA= 86° Andrews et al. 2009). The 
X 2 probability distribution is sampled by varying the 
free parameters, generating a model state and compar- 
ing this state to the data in Fourier space, following the 
Markov Chain Monte Carlo (MCMC) method described 
in Isella et al. (2009). 

For this initial model fit, we select the maximum grain 
size (a max ) and particle-size distribution slope (q) that 
best fit the unresolved SED from sub-mm to cm wave- 
lengths. For AS 209, these correspond to q = 3.5, 
&max — 1-3 mm, resulting in /3 = 1.0 between 0.88 
and 8.0 mm. Then, assuming a constant dust opacity 
throughout the disk, each wavelength is fitted separately 
with the aforementioned physical model. We obtain best- 
fit S(i2) and T(R) profiles, that will be different for each 
wavelength if the assumption of a radially constant dust 
opacity is not valid (Isella et al. 2010). 

Figure 2 shows the best-fit and 3<r constraints on T(R) 
and the disk optical depth r\(R) = k,\E(R), inferred 
from modeling each observation separately. As can be 
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Fig. 3. — Continuum emission towards AS 209, observed at 0.88 mm (top panels), 2.8 mm (middle panels) and 8.0 and 9.8 mm (combined 
through multi-frequency synthesis, bottom panels). Each observation, accompanied by the best-fit disk emission and a residual map 
obtained by subtracting the best-fit model from the observations, used Briggs weighting with robust = 0.7 (SMA, CARMA), while VLA 
data used natural weighting. Contours start at —3a, stepping by 3a (CARMA, SMA) and 6a (VLA), where a is the RMS noise on each 
map: <tsma = 4.4 mJy/beam, ccarma = 0.47 mjy/beam, <7vla = 0.01 mjy/beam. 



seen in this figure, a differing S(i2) and T(R) was inferred 
for each wavelength when k,\ is assumed to be constant 
with radius. We constrain R t to be large for the short- 
wavelength emission (R t = 61 ± 2 AU from 0.88 mm 
and R t = 60 ± 6 AU from 2.8 mm), compared to R t 
inferred at long-wavelengths (R t = 24 ± 3 AU from 8 
mm, R t = 26 ± 4 AU from 9.8 mm). We note that the 
derived temperature profiles differ by < 4 K for R < 
70 AU, and by < 10 K for R > 70 AU. The best-fit 
model visibility profiles at each wavelength are shown 
as continuous lines in Figure 1. Maps of the observed 
emission, with corresponding best-fit model and residual 
emission maps are shown in Figure 3. 

Since the disk is optically thin (right panel, Figure 2), 
the observed emission will depend directly upon the dust 
opacity and disk mass: S\ cx ft a x £ x B\(T). Hence, 
our observations constrain the product ^xSx B\(T), 
where all these physical parameters (k^,E,T) depend 
on the radius, but only the dust opacity depends on the 
observed wavelength. The modeling presented above, as- 
suming a constant dust opacity with radius and fitting 
each wavelength separately, has resulted in a wavelength- 
dependent T>\(R) and T\(R), which is unphysical but 
telling: the assumption of a constant K\ with radius is 
likely incorrect. Radial variations of the dust opacity 



are required to reconcile the differing visibility profiles 
(Real part vs. distance, Figure 1) and the wavelength- 
dependent E A (ii) and T X (R) (Figure 2). 

5. DISCUSSION 

5.1. Radial variations of f3 

Expanding on the work presented in Isella et al. 
(2010), if the true but unknown disk physical quantities 
are ~fi\(R), S(i2), and T(R), then at each wavelength 

ft a E A (i2) Bx(T x (R)) = kx(R) B X (T(R)) (2) 

where the right-hand side of this equation corresponds to 
the true disk properties, and the left-hand side encom- 
passes the constraints found from our multi-wavelength 
observations. Taking the ratio of equation 2 for Ai 7^ A2, 
and assuming that at long wavelengths k,\ oc A _/3 , we 
obtain 



Ap(R) = log" 1 



"Ai" 


x log 









E Al B Al (T Al )/B Al (T) 



^a 2 B\ 2 (T\ 2 )/ Bx 2 (T) _ 



(3) 



where Sa, T\, and T depend on the radius R. Equa- 
tion 3 is a useful prescription to constrain radial varia- 
tions of the opacity spectral slope (Af3(R) = f3(R) — f3c) 
in dual- wavelength observations. Here /3c corresponds to 
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Fig. 4. — Left: Dust opacity spectral slope, /3, vs. radius, inferred from mult i- wavelength observations of the AS 209 disk. Black line: 
best-fit /3(R), colored areas: confidence interval constrained by our observations. Vertical dashed-lines indicate the spatial resolution of our 
observations, errorbar in top-left corner indicates additional systematic uncertainty on (3(R) arising from amplitude calibration uncertainty. 
Right: Dust opacity (normalized at 300 GHz) for a m ax between 0.1-10 cm. Note that the power-law assumption, k u oc v@ , breaks down 
for (sub-)mm-sized grains. 



the spectral slope of the assumed radially constant dust 
opacity. 

In the Rayleigh- Jeans regime (when hc/X <C 
ks' Boltzmann constant), knowledge of T(R) to derive 
A/3(R) is unnecessary since B X (T) oc T. However, in the 
cold outer disk and at short wavelengths, the Rayleigh- 
Jeans assumption is ill-founded (e.g. hv ~ for 
T = 20 K at 0.88 mm). Consequently, we require an 
estimate of T(R). Since the temperature profiles inferred 
from each wavelength (T\(R)) are not very different, we 
average them to deduce T(R) (dashed- line on Figure 2). 

Equation 3 shows that in logarithmic space A/3(R) is 
the slope of a line that goes through points {x = log(A), 

y = log[ExB x (T x )/Bx(T)]}, with y evaluated at ra- 
dius R. Hence, to measure Af3(R) we find the slope 
of a straight line - in a least-square sense - that fits our 
constraints on the product Ea x B x (T x )/B x (T) from 
our set of mult i- wavelength observations. To set confi- 
dence intervals on A/3(R), we make use of the Bayesian 
approach upon which our modeling is based. Since the 
result of the MCMC algorithm is a fully sampled poste- 
rior probability distribution function (PDF) for all the 
model parameters, we construct a PDF of the product 
Ea x B\(T\)/ B\(T) at each radius R and for each wave- 
length. Random samples of these PDFs are taken at 
each wavelength, the slope of the line through points 
{x = log(A), y = log[E A x B X (T X )/B X (T)}} at a radius 
R, is computed. This slope corresponds to one random 
sampling of Af3(R). Hence, the PDF for A/3(R) is con- 
structed by performing a large number of random sam- 
ples. The peak of this PDF is the best-fit value of A/3 
at radius R. Confidence intervals are derived from the 
region of the PDF that contains 68.3%, 95.5% and 99.7% 
of all samples at equal probability (la, 2a, and 3a). 

Figure 4 presents the constraints on the radial varia- 
tion of j3 obtained from our multi-wavelength observa- 
tions. The values of f3 allowed by our observations are 
significantly different than /?jgjy[ ~ 1.7, for R < 70 AU. 
Furthermore, we find a gradient on f3(R) inconsistent 



with a constant value at the 10a level. 

5.2. Radial variations of amax 

To derive equation 3, the assumption of ft a oc A _/5 
must be satisfied. We caution that for a max ~ 0.1-1 mm 
this approximation may break down (Draine 2006), as 
illustrated in Figure 4 (right). Therefore, rather than in- 
ferring a rnax (R) from f3(R), we constrain it directly by 
fitting a specific dust opacity ~fi x to the constraints on 
the product ft a x Ea x B x (T x ) at each radius. With a 
knowledge of T(R), and for a fixed set of dust proper- 
ties (composition and grain-size distribution), we esti- 
mate the values of am ax and E that satisfy equation 2, 
now written as 



S A (ii) 



B X (T X (R)) ^ 



Bx(T(R)) 



(R)) S(i?) (4) 



where the right-hand side corresponds to our model (with 
parameters a max and E), and the left-hand side has been 
constrained by our multi-wavelength observations (i.e. 

we have a PDF for the product ft A x E A (i?) x ^ a( ^ )) ). 

At each radius R, we constructed a two dimensional 
grid of parameters {a ma;r ,E}. At each point in the 
grid, we compute the product K X (a max ) x E at each 
wavelength. We then find the probability that such 
measurement will have occurred (given our observa- 
tional constraints on the left-hand side of equation 4), 
and construct the likelihood function of the parameters 
{umax,^}- Best-fit values for a max and E are found by 
maximizing the likelihood, confidence intervals are ob- 
tained from the marginalized likelihoods. 

Figure 5 presents our constraints on a max (R) and E(i?) 
for two representative values of q. The same compo- 
sition presented in §4 is assumed, however a different 
dust mixture will influence the derived am ax (R). For a 
composition that includes updated oxygen abundances 
(Asplund et al. 2009), a smaller a max (R) is inferred (x2 
smaller), well within the uncertainties of this derivation 
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Fig. 5. — Surface density (top) and maximum grain size as a function of radius (bottom), for two grain size distributions: q = 3.5 (left), 
q = 3.0 (right). Black line: best-fit, shaded region: 3a confidence interval. Assumed grain properties are specified in figure legend. We 
compare our observational constraints with theoretical grain evolution models (Birnstiel et al. 2012), that include fragmentation and radial 
drift. 



given other unknown parameters (e.g. particle-size dis- 
tribution slope). Across the disk, grains have grown at 
least up to ~ 0.5 mm, with small grains present in the 
outer disk and large grains in the inner disk. 

We compare our observational constraints with theo- 
retical models of grain growth, employing the approxi- 
mations presented in Birnstiel et al. (2012) for the evo- 
lution of a max with radius. Since the true pro- 
file constrained by our observations depends on the as- 
sumed value of q (see Figure 5), these theoretical pre- 
scriptions will depend on q as well. Our observational 
constraints on a ma x are consistent with a radial drift lim- 
ited grain population, where the head-wind felt by dust 
particles makes them spiral towards the star. A frag- 
mentation dominated population, where the turbulent 
relative motion of particles causes collisions that either 
grow or fragment these grains, seems incompatible for 
standard values of the turbulence parameter (at = 0.01, 
Shakura & Sunyaev 1973), fragmentation threshold ve- 
locity (u t = 10 m/s, Blum & Wurm 2008), and 100:1 gas- 
to-dust ratio. These curves represent barriers that pre- 
vent further size increase, hence, the smaller of the two 
is considered the upper-limit to growth. However, the 
parameters that go into deriving the fragmentation limit 
are very uncertain, while the physics and parameters in 
radial drift are better established. Either no fragmenta- 
tion barrier exists (for a low-turbulence disk, a t < 0.001, 
the maximum collision speed never reaches the fragmen- 



tation threshold velocity in AS 209, making fragmenta- 
tion impossible) or the gas-to-dust ratio is ~ x 10 larger, 
allowing for the fragmentation-limited barrier to be con- 
sistent with our a max (R) constraint. It remains to be 
demonstrated whether these models reproduce the ob- 
served P(R) profile and the millimeter fluxes. Such anal- 
ysis requires numerical simulations starting from the disk 
parameters derived here and will be discussed in a future 
paper. 

5.3. Summary 

We have obtained multi-wavelength spatially-resolved 
observations at 0.88, 2.8, 8.0, and 9.8 mm of the 
AS 209 circumstellar disk. These observations reveal 
a wavelength-dependent structure, explained as radial 
variations of the dust opacity across the disk. We find a 
change in A/3 > 1 between the inner (~ 20 AU) and outer 
(~ 120 AU) disk, inconsistent with a constant f3 value. 
This gradient in f3(R) implies that a significant change 
in the dust properties as a function of radius must exist. 
We interpret this gradient as a decrease in the maximum 
grain size with radius: going from several centimeters 
or more in the inner disk to sub-mm-sized grains in the 
outer disk. When compared with theoretical models of 
dust size evolution (Birnstiel et al. 2012), we find that 
our observational constraint on a max (R) agrees with a 
radial-drift-dominated population, for reasonable values 
of the composition, particle-size distribution, and disk 
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properties. 
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